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ABSTRACT 


An analytical study was undertaken to determine two- 
dimensional wall conduction effects during film condensation 
on an internally finned, rotating heat pipe. An earlier 
Finite Element method was modified to solve the problem 
by uSing different numbers of axial increments and different 
numbers of elements. Resulting heat transfer rates were 
Shown to be Similar to those obtained earlier. 

A rotating heat pipe was tested uSing various copper 
condensers including two smooth-wall cylinders anda 0.5 
Gegqree Eruncated cone. All comdensers were tested at 
different rotational speeds using distilled water as the 
working fluid. In the case of the 1.46 inch diameter 
cylindrical condenser, data were also taken with Freon 113 
and ethanol. The heat transfer rate of each condenser was 
plotted against the saturation temperature of the vapor. 

Using the experimental results obtained in this thesis, 
together with the data of earlier theses, a comparison 
was made with the theoretical analysis of Roetzel. Agree- 
ment is reasonable, and Roetzel's analysis can be used to 
Predict Overall rematwmg Heat pipe performance to within 
+ 20 percent for water. The data for Freon 113 and ethanol, 
however, fall approximately 35 percent higher than this 


theory. 
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io LE RODUCTION 


Ae ) Poe ROTATING HEAT PIPE 

The rotating heat pipe is a closéd rotating container 
capable of transferring large amounts of heat from rotating 
machinery. It consists of three main parts: (1) a cylin- 
@rical evaporator, (2) a cylindrical or truncated-cone 
condenser and (3) a working fluid. The operation of a 
rotating heat pipe is shown schematically in Figure l. 

When the rotating heat pipe is rotated above a critical 
speed, the working fluid forms an annulus in the evaporator. 
Heat input to the evaporator causes some of the working 
fluid to vaporize, and this vapor flows into the condenser 
transporting latent heat of vaporization. External cooling 
of the condenser allows the vapor to condense on the inner 
condenser wall. The centrifugal force due to the rotation 
has a component acting along the condenser wall which 


accelerates the condensate back to the evaporator. 


See curr iNG HEAT PIPE LIMITATIONS 

Ballback [1] has studied the limitations imposed on 
the rotating heat pipe by various fluid dynamic mechanisms 
in a preliminary theoretical analysis of its operation. 
By using existing theories and experimental correlations, 
he was able to estimate the sonic limit, the boiling limit, 


the entrainment limit and the condensing limit of performance. 


18 
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A typical rotating heat pipe with the physical dimen- 
sions as shown in Table 1 was chosen in order to demon- 


Strate the predicted performance. 


Table l 


Specifications of a Typical Rotating Heat Pipe 


Length 14.000 inches 
Minimum diameter 2720.00 inches 
Wall thickness Diel25 inches 
Internal half angle 1.000 degree 
Rotating speed 2700 RPM 


ieviesomic Lami 

As the heat flux in the rotating heat pipe is 
increased, it 1S possible that the flow rate of the vapor 
down the container may be limited by the attainment of a 
choked flow condition for the vapor. This choking of the 
vapor flow will definitely limit the amount of energy the 
vapor can transport down the pipe, thus limiting the opera- 
tion of the pipe. The heat transfer mechanism employed in 
the rotating heat pipe depends on the latent energy of the 
working fluid. The heat transfer rate in the pipe may 


then be expressed as: 
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where 
m = mass flow rate of the vapor, lbm/hr. 


V 


From continuity, the mass flow rate of the vapor may be 


expressed by: 


n =O uA (2) 


where 


& 
il 


velocity of the vapor, ft/sec 


A = cross section area for the flow, Ft’. 
Then 
Cee = NEO. >A h , (3) 


The vapor velocity is taken as the sonic velocity, so 


eMac: 
Ue Vg, KRT (4) 
where 
Cmemecon@e Velocity, £t/sec 
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Ee=1om 


ao, = geeeteration of gravity, 32.1739 IEG 
k = the ratio of specific heats, 

R = gas constant, ft-lbf/lbm-deg R 

pi = absolute temperature, deg R 
Using equations (3) and (4), the sonic limit of a 


typical rotation heat pipe can be computed, as shown in 
Figure 2. 
ee Bolling Limit 
Kutateladze [2] has postulated that the transition 
from nucleate boiling to film boiling is purely a hydro- 
dynamic phenomenon. Based on this assumption he has 


obtained the correlation formula for predicting the burnout 


heat flux: 

; 1/4 

Sto S/S heg {oglo~ - 9,1) (5) 
where 

AL = heat transfer area in the boiler, Ft 

0, = density of the vapor, Lbm/£t> 

ae = latent heat of vaporization, BTU/lbm 

fe} = surface tension, lbf/ft - 


ae, 
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Figure 2. Operating Limits of a Typical, Water-Filled, 
Rotating Heat Pipe 
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= acceleration of gravity, ft/hr- 


density of fluid, Usa Be> 


constant value. 


Kutateladze's experimental data indicate a value of 


K between 0. 


tse. Using equation (5), the boiling 


limit of a typical rotating heat pipe can be computed, as 


shown in Figure 2. 


3. Ehemsanmene, Gbimit 


Sakhuja [3] studied the flooding constraint in 


wickless heat pipes, where he assumed that: 


a) 


b) 


Se) 


d) 


e) 


the vapor and liquid are in counterflow motion 
and are exposed to each other; 

Ee reingertim Of Ligurd (under gravity) is 
subjected to a shear stress by counterflowing 
vapor, which tends to retard the film; 

the liquid film remains smooth and stable, and 
the shear stress is usually small; 

the buoyancy forces, due to density differences 
between the vapor and the liquid, are responsible 
for maintaining the counterflow; and 

the buoyancy forces are balanced with dissipative 
effects which are proportional to momentum 


Fluxes of vapor and liquid streams. 


He has obtained a correlation for flooding: 
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where 


qd, = heat transfer rate, BTU/hr 
A = flow area, aya 
x 
C = dimensionless constant, 0.725 for tube 


with sharp edged flange 


D = inside diameter of the heat pipe, ft 
g = acceleration due to gravity, ft/hr- 
o. = density of the fluid, Lbm/ft> 

e.. = density of the vapor, lbm/£t> 

as = latent heat of vaporization, BTU/1lbm. 


Using equation (6), the entrainment limit of a 
typical rotating heat pipe can be computed and is shown 
mee rigure 2. 
fe COndensingmolimat 
The solution for condensation in the rotating heat 
pipe has been carried out by Ballback [1]. He modeled the 
condenser section of the heat pipe as a rotating, truncated 


cone, and an equation to compute the condensing limit was 


derived: 
. Ke oe w he (Tr -7 > 27" 8/3 8/3. 3/4 
qd. = nf2 £0E* Neh sw {{R, + Lsing] - RS J (7) 
u sin¢ 
fe 
where 
qd, = total heat transfer rate, BTU/hr 
Ke = thermal conductivity of the condensate 


fon, eolU/Ne=rt-daeq F 
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= density of fluid, lbm/ft? 


Pe 

W = angular velocity, l1/hr 

Bae = latent heat of varporization, BTU/lbm 

T = saturation temperature, deg F 

T. = inside wall temperature, deg F 

He = viscosity of fluid, lbm/ft-hr 

) = half cone angle, degrees 

Ro = minimum wall radius in the condenser 
section, ft 

L = length along the wall of the condenser, ft. 


Equation (7) is a function of the physical proper- 
ties of the working fluid, the geometry of the rotating 
heat pipe and the operating speed. Using this equation, 
the condensing limit for a typical rotating heat pipe can 
be computed. This result is also shown in Figure 2. 

It is clear from the results of Figure 2, that for 
the geometry and working fluid chosen, and under normal 
circumstances, the heat transfer rate is limited primarily 
by the amount of heat which can be transferred through the 
condenser wall. However, other limitations may become 
important as heat pipe geometry and operating conditions 
are altered. 
©. §HEORETICAL INVESTIGATIONS OF FILM CONDENSATION 

ON SMOOTH ROTATING CONDENSER SURFACES 
Sparrow and Gregg [4] were the first to study condensa- 


tion on a rotating surface. They studied condensation on a 
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rotating disk at uniform surface temperature, se Situated 
in a large quiescent body of pure saturated vapor. They 
assumed that the condensed liquid formed a continuous film 
on the disk and that fluid in this film moved radially 
outward along the disk under the action of the centrifugal 
force field. They arrived at the following solution for 

the situation where energy convection and acceleration terms 


are negligible: 


v 2 
h (2) */ Pr 
eee Ot ee (8) 
i p fg 
where 

h = local heat transfer coefficient, 
BTU/hr-ft2-deg F 

K ¢ = thermal conductivity of condensate, 
BTU/hr-ft-deg F 

v = kinematic viscosity of condensate, GF Jive 

W = angular velocity, l/hr 

Pr = Prandtl number, ee 

u = absolute viscosity of condensate, lbm/ft-hr 

Ce = specific heat of condensate, BTU/lbm-deg F 

aie = disk surface temperature, deg F 

ae = saturation temperature of vapor, deg F 

avlk = ar = Toe deg F 

ae = latent heat of vaporization, BTU/I1bm. 


Ze 





Equation (8) predicts that the local heat transfer 
coefficient is constant across the disk and does not vary 
Wie rad@al position. 

Sparrow and Hartnett [5] continued the research 
using a rotating cone with the same assumptions used earlier. 
Their result for the heat transfer coefficient was: 


NA 
(sin 6) (9) 


Dooncl Sagar 
where ¢ is the internal half-angle of the cone. As pointed 
out by the authors, however, equation (9) applies only to 
cones that are not too slender. Using their result, the 


mean Nusselt number for a cone can be written as: 


Nu eo ae 0 ome ace Le Pr eae 
m k ° Oe yn 
Vv p fg 
(250) 


Leppert and Nimmo [6,7] studied laminar film 
condensation on a surface normal to body or inertial forces 
Such as a finite horizontal surface, or a rotating drum. 

They found that the local heat transfer coefficient varies 
with position, and the mean Nusselt number, for a constant 
Surface temperature, with negligible surface tension effects, 
1S% 


IAS 


Nu = 0.64 Sh (i) 
™m 


Ze 





where 





le GL 
Nu, = mean Nusselt number, a 
g 97 Deg ie 
Sigh = Sherwood number, Soa 


Equation (11) can be applied to the case of 
condensation on the inside of a rotating cylinder of radius, 
r, where g is replaced by Puen 

The first theoretical investigation conducted at 
the Naval Postgraduate School was performed by Ballback 
[1]. His analysis was a Nusselt type analysis for film 
condensation on the inside surface of a rotating truncated 
cone of half angle cone, $¢, with a geometry defined in 
Figure 3. His assumptions were Similar to those used in 
the Nusselt analysis of film condensation on a vertical 
wall, namely: 

1) steady state condition, 

2) film condensation, not dropwise condensation, 

3) laminar flow of condensate, 

4) constant fluid properties, 

5) no condensate subcooling, 

6) no momentum change through the condensate, 

7) no liquid-vapor interfacial shear stresses, 


8) linear distribution of temperature through the 
condensate film, 


9) negligible force of gravity, and 


10) constant wall temperature. 
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Additionally Ballback assumed: 
11) avery small film thickness, 4(x), 


12) no radial distribution of pressure in the 
vapor space, 


13) negligible density of the vapor, 


14) negligible velocity gradients in the 
circumferential direction, and 


15) very small slope of the condensate film, 
dé/dx when compared with tan $6. 


He used the boundary condition, that 6(x) must be 
zero at x = 0, and arrived at the following solution for 


the local film thickness: 


D3 ke Ug(T, - T,) ae ees 4 — 
2 con. R,+x sin ¢ 
f ney 
where 
6(x) = film thickness, ft 
Ke = themmal cConauctivity of the fluid, 
Buu7he-fE-deg F 
Ue =Evyiscosity Of the fluid, lbm/ft-hr 
o, = density of fluid, lbm/ft* 
Ww) = angular velocity, l/hr 
ae = inside wall temperature, deg F 
Te = saturation temperature, deg F 
eg = latent heat of vaporization, BTU/1bm 
d = half cone angle, degrees 
Ro = minimum wall radius in the condenser 


section, ft. 
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The local heat transfer coefficient, h(x), is 


given by: 


Baec- 2 
asin Ae 1 


I ————————— 
3S Ak T a4 Ro ae 


Roo x Sim 5! 
O 


1/4 


(13) 


The mean heat transfer coefficient over the whole 


surface area of the truncated cone is given by: 


L 
: 1 | 
ha Shae 2R_ +Lsin 3 Pd 2mh(x) [R +x San sodieodx 


(14) 
where 
L = length along the wall of the condenser, ft. 


Integrating equation (14) yields a value for the 


mean Nusselt number: 





a L ae ie Re oe 1/4 
Nu = = 0.904 (——~— = ) G(B) 
I< y2 Co be qT, fg 
(5) 
where 
8/3 a 3/4 


78 (2+ 8) 


Zi 





and 


Tr sin > 


Bemecron (15) can be put into the form: 


Nu, = (NU) on —— CL), 


where (Nu) og is the result of Sparrow and Hartnett [5] 
for a cone, as given by equation (10). The above result 
was confirmed as a special case of the later work of Dhir 
and Lienhard [8], who studied laminar film condensation on 
axisymmetric bodies in nonuniform gravity fields. 

Daley [9] was able to improve on the theory intro- 
duced by Ballback by consideration of the thermal resis- 
tances in the condenser wall and on the outside surface. 
Both Ballback and Daley, however, neglected the drag due 
to vapor friction on the counterflowing condensate film. 

Newton [10] extended the analysis performed by 
Ballback and Daley to include the effects of vapor pressure 
drop and interfacial shear stress between vapor and conden- 
Sate. Using the same geometry in Figure 3, he assumed that 
the momentum changes within the condensate film were small 
and arrived at the following differential equations for 
condensate film thickness, vapor velocity and heat transfer 


rate: 
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2 P, (x) Po (x) WER (x) ae 
Py. PS EN (x) i a a, ae OGY (x) ——.o 2P, (x) ao] 
Ss = tant) > 
dx 

5 20 v(x) 
P, (x) cos¢ [p .w R(x) + ae 7 
2(R_+xsin 9d (T ad ee 
oe SC} NBO tnd > Gee, 0) 
ax 2 R(x) ax 
(eee hye. R(x) 
V Ke k gle ike 
(19) 
a ei aa (20) 
ax 6 (x) in ty AL 
Ki k h 
1 W 
where 
Baise) 5 4 
P, (x) = (R, + xsin $ ls - 57C0s OMOnS) Wis, 
2 
: 6 (x) 
P. (x) = (R_+xsin >? )—=—— - cos 6(x) , 
Oo Zz ene 
v (x) = velocity of vapor, ft/hr 
Ae = saturation temperature of vapor, deg F 
Te = ambient coolant temperature, deg F 
2 = thickness of the condenser wall, ft 
Bat) = film thickness of the condensate, ft 
Kp = thermal conductivity of condensate film, 


BTU/hr-ft-deg F 
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k = thermal conductivity of condenser wall, 


Y BTU/hr-ft-deg F 
h = external heat transfer coefficient, 
BTU/hr-ft*-deg F 
Fag = latent heat of vaporization, BTU/1lbm 
£ = friction factor, dimensionless, 
R, = vapor Reynolds number, 2R(x) P v(x) /H 
R(x) = radius along the condenser wall, 
Ro + xSingd- dcos 
) = half cone angle, degrees 
, = density of vapor, ery ee 
Le = viscosity of vapor, lbm/ft-hr. 


The integrations were started by assuming initial 


conditions 6 = Ft and v = 0 at x = 0. With these initial 
eendaitions s = tangd at x = 0. Newton obtained results 


per nalt cone angles of 0, 0.1, 0.2 and 0.3 degrees. As 
the half cone angle increased, the solution of the differ- 
ential equations became extremely sensitive to the initial 
values of 55. Hewalseenaceauthnat his results ford (x) 
Showed that (except near the ends of the condenser) 

¢ << tan >. He, therefore, made this approximation and 
Simplified the model further using an order of magnitude 
analysis of the terms in equations (18) and (19). With 


these simplifications, equations (18) and (19) were 


reduced to: 
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ay) _ 2(R. + xsin aT. =i 


2v(x)sin 6 





ax — oo (21.3 

re a a R(x) 

‘ox. if K + A) £5 R(x) 

2 
st : 20 v(x) 2 on 
———R (x) sin o[0 0 15 5.4) as =S! = pv (x) ETR(x) r ] 
2 
se) 


Equation (21) was numerically integrated using the 
initial condition that 55 = 0 and v = 0 at x = O and using 
miemecuite eCuatron for 6(x), equation (22). After the 
integration was completed for §(x) and v(x), the energy 
equation, equation (20), was integrated over the length of 
the condenser to find the total heat transfer rate of the 
heat pipe. This theory was compared to the experimental 
results of Schafer [15] and Tucker [20], predicting perfor- 
mance to within +20 percent for water, ethyl alcohol and 
Peaeon 13. 

Daniels and Al-Jumaily [11,12] extended the analysis 
Of Ballback by including drag effects due to the vapor 
fewecion. Their solution for the local film thickness 


distribution along the condenser is given by: 
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where 


a 


0* (wr-g) sinoh, x? 
Shy = Sherwood number, TuKAT 
2 
. OT vo Eg* cos 
Dry. = drag friction number, OKAT 
pu. x cos 6 
Re = two phase Reynolds number, ae 
Vx u 
0 = liquid density, lbm/ft> 
piowe= es Elim thickness, £t 
W = angular velocity, l/hr 
d = half cone angle, degrees 
T = shear StKESS at liquid vapor interface, 
sf Wess etm 
u = liquid condensate velocity, ft/hr 
g = acceleration of gravity, ft/hr- 
a = latent heat of vaporization, BTU/1bm. 


At any section x in the condenser, since heat is 
transferred through the condensate layer by conduction 
only, the local heat flux Q. = kAT/6 (x) and equation (23) 
can be re-arranged as follows: 

i a ul sas 3 kAT, 4 


4 
Q. - ea ee Qo. + 3Dr (——) Q.. - Sh, SS =P (24) 


a2 





where 


Q, = local heat flux, BTU/hr-£t*. 


This equation relates the three main parameters 
controlling the rotating heat pipe operation assuming 
fixed geometry and working fluid. The mean Nusselt number, 


Nua for this analysis was, 


where suffix L refers to quantities evaluated at the end 
of the condenser section x = L. 

Marto [13] compared the heat transfer rates of 
Daniels and Al-Jumaily with Ballback and the approximate 
numerical aiken £o equataoms (20), (21) amd (22). For 
the case of Freon 113, equations (20), (21) and (22), 
which include both effects of interfacial shear as well as 
vapor pressure drop, gave slightly lower heat transfer 
rates than the analysis of Daniels and Al-Jumaily which 
includes only interfacial effects. Also the analysis of 
Bupeedeckepredicts slightly larger results than the computer 
analysis but lower results than those of Daniels and 
Al-Jumaily. This may be due to differences in the Freon 113 
fluid properties used or due to the geometrical inconsistency 
in matching all the geometrical dimensions as stated in the 


analysis of Daniels and Al-Jumaily. For the case of water, 
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Ballback's results again predict a slightly higher heat 
transfer rate than the approximate numerical solution, due 
to his neglecting effects of interfacial shear and vapor 
Pmesot~e creep. When these effects are neglected in the 
approximate numerical solution, the results agree very 
well with those of Ballback. For these reasons, the 
effects of interfacial shear and vapor pressure drop can 
safely be neglected. 

Roetzel [14] studied condensation in steam heated 
fast rotating paper drying drums, as shown in Figure 4. 
In his analysis, the heat transfer is determined mainly by 
the method of condensate removal. The condensate is collected 
in a circumferential groove on the inside surface, and from 
there it is removed by means of a stationary syphon with 
its tip in the groove. The mean film thickness is slightly 
Smaller than with a rotating syphon because the condensate 
is sucked off from a lower level. The solution for the 
case of a drum with an infinite radius of curvature (The 
inside surface is exactly cylindrical.) with a uniform 
temperature difference is given by: 


75 
Nu, = (0.86 E) (26) 


where 


Nu. = Nusselt number for the cylindrical drum, 
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Figure 4. Schematic Drawing of Paper Drying Drum 


with Two Curved Conical Sections and 
One Sampling Groove 
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ben, L 
fee OGK AT 
ieeern 
Ke = thermal conductivity of condensate, 


Bed/7ar-ft-deg F 


Sh = Sherwood number 

oo = acceleration normal to the wall, ft /hr- 

O = density of condensate, Lbm/ft> 

tag = latent heat of vaporization, BTU/lbm 

L = total length of condensate, ft 

Vv = kinematic viscosity of condensate, Bee Pee 


Diticwecqtattion can be written in terms of the 
Sherwood number: 


1/5 (27) 


Nu_ = Oey om Sh) 
and compares favorably to the analysis of Nimmo and Leppert, 
equation (ll). 

For the case of finite radius of curvature, Roetzel's 
solution becomes: 


OLE 2277156 


Nu. = Ga) + (0.86E) 


237 2074723 
oF ] 


(28) 
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where 


P= rte tuS GO. Curvature, ft, and 


ot 
ll 


E@Qual flow=lengehn of condenser, ft. 


ties eGUuatiton Can be put in terms of Sh, Nu,, R, and L as: 


L) 23/16 


it 23/80) 4/23 
R 


Nu = Nu, fl +1.924 ( Sh 28h) 
Thus the Nusselt number can be directly determined when 

the ratio L/R and the value of the dimensionless group, E, 
are known. For the truncated cone, rotating heat pipe, as 
shown in Figure 3, equation (29) can be applied to compute 


the heat transfer by defining the radius of curvature 


R= L/tand; the solution then becomes: 


Nu = Nu, (1 +1.924 tang °7/4° gn*9/8914/23 (39) 
where 

W) = half cone angle, degrees 

Nu. = Nusselt number of the cylindrical drun, 


equation (ll). 


Equation (30) has never been compared with experi- 
mental data, so using the existing rotating heat pipe data 
together with new experimental data at different half cone 


angles will provide a check on its usefulness. 
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These theoretical investigations are summarized 

below: 

1. The vapor pressure drop and interfacial shear 
stress effects between vapor and condensate can 
be neglected. | 

2. The analysis of Roetzel [14] appears to be the 
most appropriate for use in a rotating heat pipe. 

3. The important parameters are: 

(a) speed, 
(b) geometry, and 
(c) Sherwood number, which includes the properties 


of the working fluid. 


D. ANALYSIS OF FINNED ROTATING CONDENSER SURFACES 

The use of axial fins on the inner surface of the con- 
denser can improve heat transfer performance through two 
mechanisms: a decrease in condensate film thickness near 
the apex of each fin and an increaSe in inner wall surface 
area. Schafer [15] developed an analytical model of the 
internally finned condenser for the case of a fin of tri- 
angular profile as shown in Figure 5, with a specific 
geometry as shown in Figure 6. In the analysis, he assumed 
one-dimensional heat conduction from the inside to the 
outside surface and a linearly decreasing surface tempera- 
ture profile from the fin apex to the midpoint of the 


trough. Results from his analysis indicated significantly 
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Figure 5. Internally Finned Condenser Geometry 
Showing Fins, Troughs and Lines of Symmetry 
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Figure 6. Condenser Geometry Considered in Analysis 
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improved heat transfer rates for the finned condenser 
compared to the smooth condenser. 

Corley [16] developed a numerical solution for this 
case using the Finite Element method to account for two- 
dimensional conduction effects within the wall for the 
triangular fin profile and a parabolic wall surface tempera- 
ture distribution along the fins. His results indicated 
Significantly improved heat transfer performance with heat 
transfer rates about 75 percent above those predicted by 
the one-dimensional linear temperature distribution analy- 
Sis of Schafer [15]. However, Corley noted that an error 
of 50 percent in the calculation of apex heat transfer rate 
might then give a total heat transfer error of as high as 
15 percent. Two possible solutions to the boundary con- 
ert1on Menien were recommended. The first was to increase 
the number of Finite Elements in the analysis, and the 
second solution was to consider a wall geometry as shown 


mer lgunue 7’. 


Peelaeots OBJECTIVES 
The objectives of this thesis were therefore: 
iy to find out which existing theory best correlates 
Rotating Heat Pipe data by comparing these theories 
with existing data and new experimental data of 
this thesis, and 
2) to further study finned condensers by improving 


the Finite Element results obtained earlier. 


4l 





Figure 7. Alternate Condenser Wall Geometry 
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Piet ite Been yT ANALYSIS AND RESULTS 


A. REVIEW OF EARLIER ANALYSES 

Schafer [15] studied the one-dimensional heat transfer 
solution, and Corley [16] studied the two-dimensional heat 
transfer solution to finned condensers. Both used the 
Same asSumptions and boundary conditions, based upon the 
work of Ballback [1], which are similar to those used in 
the Nusselt analysis of film condensation on a vertical 
wall. The more important of these assumptions are listed 
below: 

1) steady state operation, 

2) film condensation, as opposed to dropwise condensation, 


3) laminar condensate flow along both the fin and 
the trough, 


4) static balance of forces within the condensate, 


5) one-dimensional conduction heat transfer through 
the film thickness (no convective heat transer), 


6) no liquid-vapor interfacial shear stresses, 
7) no condensate subcooling, 


8) zero heat flux boundary conditions on both sides 
of the wall section, as shown in Figure 6, 


9) saturation temperature at the fin apex, 
PO) zero riim thickness at the fin apex, 

ll) negligible curvature of the wall. 

The complete development of Schafer's finned condenser 


heat transfer theory is presented in his thesis [15]. His 
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results are used for comparison purposes in this 
thesis. 

Corley's study [16] considered an internal fin geometry 
as shown in Figure 6. Two-dimensional conduction through 
the condenser wall was considered in the heat transfer 
analysis by using a Finite Element solution. Instead of 
using a linearly decreasing temperature profile along the 
wall from the apex of the fin to the midpoint of the 
channel as Schafer did, a parabolic distribution was assumed 
to calculate average nodal convective coefficients. The 
eight node finite element mesh used in the two-dimensional 
analysis of the axially finned condenser is shown in Figure 
8. A computer program for the case of two-dimensional con- 
duction was developed using a modified version of the three 
dimensional, isoparametric Finite Element program assembled 
by Lew [17]. 

A parabolic temperature distribution is assumed to exist 


along the fin surface using the boundary condition, T= T 


S 
at z = 0, since 6(0) = 0. Therefore; 
T(z) = a.z~+b.24+T 
W 1 i S 
or 
TM (2) = -a.z- - b.z (31) 
S W l Jk 
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Figure 8. Condenser Geometry Considered in 
Two-Dimensional Analysis with 2 Finite 
Elements 
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where 


ay and Dy are unknown constants, and 


z is the distance along fin surface from apex of 
fone £t. 


From Schafer's analysis [15] a differential equation 
for fin condensate thickness, 6(2), can be derived: 
Ke(TL - T(z)) Lys Z 


5° (z) a6 (z) 2 ee (2) 


208 ; 
Pw (R, + xsin ee cos ¢ cos a 


where 


6(z) = the fin condensate thickness, ft 

Ke = thermal conductivity of working fluid, 
BHM) gle Seeley, 19 

Pe = density of working fluid, lbm/f£t> 

Ue — iemcececsityeOrevorcing fluid, lom/ft-hr 

W) = J sangular velocity, i/nhr 

Ry =minimnegm radius; ft 

x = distance along apex of fin, ft 

Z = distance along the fin, from the apex, ft. 

d = half cone angle, degrees 

fo} = fin half-angle, degrees 

See = latent heat of vaporization, BTU/1lbm. 
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Equation (31) is then substituted in (32), and the 


result is integrated from z = 0 to z to give an equation 
Gor 0 (2) : 
23 22 
4k.(-a, = -b, =) U 1/4 
Ga _ 3 ie f (33) 
PW (Ro + xsin 2) cos 0 cos a 


The local convective heat transfer coefficient assuming 
one-dimensional heat conduction through the thin condensate 


film is then: 


k k2 pu (R +xsin gdh cos¢cosa 1/4 
ie) _ f = a ak O fg 
fin On za) 3 2 
Aes = b SH) 
if ee ILM 
(34) 
which is applicable from z = 0 to z* = Z, 7 S*/cosa, 


where 6* is the trough condensate thickness, as shown in 
Figure 8. Across the width of trough, from z = Zo to the 


mid channel position, the heat transfer coefficient is: 


*e 


hiGeyveugh OF tx) (35) 


where 6*(x) 1S considered constant across the trough width. 
Based upon a check of the wall temperatures generated from 
Schafer's computer program, the fluid temperature difference 


across the width of the trough was predicted to be small. 
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Fluid thermal conductivity values were assumed to remain 

Relacively Comical t aeetSea result, h(x) 1s assumed 
trough 

to be constant across the trough width and is calculated 

using equation (35) and a representative average inner 

wall temperature to determine Ke. The convective coeffi- 


cient along the outside surface of the wall, h is also 


Oise 
assumed to be constant. The modified trough mass flow 


equation becomes: 


yo | a: 
OW (R, + xsin ¢)5* (x) sin db 2 


M (x) = Janne | oe TOs (x) tan o| : 
£ 


(36) 


Within the computer program, the heat transfer rate 

for the entire condenser is determined by dividing its 
axial length into one hundred equal length increments of 
length, Ax (See Figure 6.), solving for the incremental 
heat transfer rate uSing an iterative procedure, and 
Summing these incremental rates to yield a total heat 
transfer rate. A flow chart description of the computer 
program used in Corley's analysis is given in Figure 9. 
At each increment, constant values of ay and Dy are deter- 
mined using a Lagrangian interpolation fit of the tempera- 
tures, at nodes 1, 3 and 6 as calculated from the previous 
iteration. For the first increment, ay and by are calcu- 


lated from an initial guess of the temperatures at these 


nodes. Convective coefficients at several locations along 
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Figure 9. 
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Computer Program Flowchart of Two-Dimensional 
Conduction Analysis of Corley [16] 
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the fin are determined using equation (34). Average 


nodal values, h , and the average h of the 


lavg’ Sa ere 
left-hand boundary region associated with node 6, in Figure 
8, are calculated using the position dependent values of 
A(Z) es, in a Simpson's Rule approximation. In the trough, 


h (x) is calculated from equation (35) using the 6* 


trough 
value calculated from the previous increment. This value 
is used as the average convective coefficient for nodes 9 
and 12 and the right hand boundary region of node 6. The 
average convective coefficients calculated for the left- 
and right-hand boundary regions of node 6 are weighted 

by their respective region lengths. The two values are 
Summed to yield the average convective coefficient for node 


S18! The remaining boundary conditions are applied, 


6avg> 
and the incremental heat transfer rate per unit condenser 
length, Qa, and nodal temperature values are generated 

from the finite element subroutines. If the resultant 

heat transfer rate converges to within the criterion of 

0.05 percent change from the value of the previous iteration, 
ee” Py-1) 0.0005 


25 


r (a7) 


then the incremental problem is deemed solved and Q. 1s 
set equal to the final iteration, Q;. he-pon the other 


hand, the convergence criteria is not met, newly calculated 
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nodal temperature values are used to establish a new set 

of inner wall convective coefficients, and the finite element 
solution is repeated until the heat transfer rate converges 
to within the acceptable 0.05 percent tolerance. Because 

Q is initialized at zero for each increment, the iteration 


procedure is repeated at least twice for every increment. 


Incremental mass flow rate is calculated uSing: 
| Q, 
AM = 2 —— Ax ; (38) 


This value is added to the mass flow rate from previous 
increments. The subsequent interval's trough condensate 
thickness, 6*(x), is calculated using a polynomial root- 
finder subroutine to solve eee (36). This incremen- 
tation is continued until the entire length of the conden- 
ser 1s traversed. Incremental heat transfer rates are 


summed, and a total heat transfer rate becomes: 


bt 
© 


O° 
Il 


0 
tot 2 Ne : (Q;. Ox) (39) 


fe. 
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where Nein is the number of fins along the inner wall 
circumference. To start the incrementation, two initial 
assumptions are made. The initial value of 6* = 5 ¥ is 
taken from the analysis by Sparrow and Gregg of condensa- 
Euemmonwa FOtLating disk [4]. In addition, at the first 


increment, initial guesses of fin temperatures at nodes l, 


3 and 6 are made in order to begin the solution. 
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Boop tle TeATrION TO THE FINITE ELEMENT ANALYSIS 

This thesis is a continuation of Corley's [16] study 
of the addition of internal axial fins to the condenser 
surface of a rotating heat pipe. He recommended that one 
way to approximate more closely the actual heat transfer 
rates would be by increasing the number of elements for 
the same geometry. Since Corley used one hundred equal 
increments of length to compute the heat transfer rates, 
the first step in this thesis was to reduce the number of 
increments to fifty and twenty five, to minimize the com- 
puter time. Using Corley's program, at a speed of 3000 
RPM, with external convective coefficients of 1000, 1500 
and 5000 BTU/hr-ft*-deg F and fin half angles of 10, 30 
and 50 degrees, the results for water as the working 
Eiwad are shown in Tables 2, 3 and 4. It is clear that 
reducing the number of increments to twenty five had little 
effect upon the heat transfer results (approximately 3 
percent error). 

To reduce any error in the calculation of the heat 
transfer rate at the apex of the fin, the number of Finite 
Elements has to be increased between the apex of the fin 
and the trough, as shown in Figure 10. The geometry becomes 
3 elements and 18 nodes. The computer program was therefore 
modified to solve this case, uSing the same boundary condi- 
tions and assumptions. These results are shown in Tables 5 


through 10. Upon dividing element (a) Mierigire 10 ante 
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Theoretical Finned Condenser Heat Transfer 


TABLE 2 


Rate Versus Number of Axial Increments 


— ome 2. 
ee = 1000 BTU/hr-ft -deg F 
Fin Half Ne. of aa 
Angle Fins 
10 276 100 
50 
200 
30 84 100 
Lo 
200 
50 40 100 
150 
200 
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44427. 


OC. 


Zibb 6. 


eSe22. 


2670. 


25255. 


42623. 


50 


POOF 


26952. 


43415. 


Teo. 


AOS. 


42627. 


7o43. 


AB a5) e 


41641. 


At 3000 RPM 


of Axial Increments 


100 


GI22. 


26422. 


42911. 


97 S0R 


ZO 949... 


21.33% 


953.0% 


293497, 


fA 





Theoretical Finned Condenser Heat Transfer 
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Theoretical Finned Condenser Heat Transfer 


Rate Versus Number of Axial Increments 
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Figure 10. Condenser Geometry Considered in 
Two-Dimensional Analysis with 
3 Finite Elements 
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two equal parts, as shown in Figure ll, the geometry 
becomes 4 elements and 23 nodes. The computer program 
was modified again to solve this case, using the same 
boundary conditions and assumptions as before. Results 


for this case are also tabulated in Tables 5 through 10. 


fae eDisCUSsiON OF NUMERIGAi RESULTS 

Computer results for the finned copper condenser were 
generated for three outside convective coefficients, three 
fin angles and two rotation rates. The working fluid used 
throughout the analysis was water. Tables 5 through 10 
compare the results of using 2, 3 and 4 elements, with those 
Seasemarer [15]. 

Using the Finite Element solution, the heat transfer 
rates change less than 10 percent, as the number of elements 
change, so the heat transfer rates in the finned condenser 
can be approximated by 2 elements, and the error in the 
calculation of the heat transfer rates at the fin apex can 
be neglected. In comparing the smooth condenser to the 
finned condenser, the most dramatic improvement occurs 
for a = 50 degrees with 40 fins. It is clear in this 


situation that Qy eee is far greater than 


smooth 


4 Element/Asmooth’ 
The heat transfer rate was found to vary slightly with 


rotational speed in a finned condenser. This result agrees 
with the results obtained earlier by Schafer [15]. Schafer 


Suggested that the invariance in heat transfer rate with 
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Figure ll. Condenser Geometry Considered in Two- 
Dimensional Analysis with 4 Finite Elements 
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the rotation rate was due to the cancelling effects of 
decreasing film thickness and decreasing temperature 
difference across the film thickness as rotation rate is 
increased. The influence of external heat transfer 


coefficient, AS is shown in Tables 5 through 10. A 


based 


large external heat transfer coefficient (Root = 5000 
BTU/hr-£t-deg F) is seen to produce the largest heat 


transfer rates. 


65 





III. EXPERIMENTAL PROGRAM 


A. EXPERIMENTAL EQUIPMENT 
1. Description of Equipment 

The equipment used to support the experimental 
program was the same equipment used by Loynes [18] and 
Wagenseil [19]. Figures 12 and 13 depict the experimental 
apparatus. The evaporator was a copper cylinder with 
O-ring seals at each end. Heat input to the evaporator 
was provided by a resistance heating element wrapped 
around the evaporator. The power supply for the heater 
element was changed from load bank resistors to a solid 
State phase amplifier power controller to reduce the space 
and inconvenience of the load banks. The model PA-1-2490 
phase/amp is a single phase silicon controlled rectifier 
power amplifier designed to control and regulate A-C power. 
A six-volt battery was used as a power supply for the 
Gemctrolier circuit. 

The condenser was surrounded by a sealed insulated 
box. Garlock seals at each end of the condenser kept the 
cooling water from leaking out. The cooling water system, 
as shown schematically in Figure 14, was four supply tubes, 
0.25 inch diameter with 0.031 inch holes drilled every 
0.20 inch. The supply tubes ran the length of the conden- 
ser, parallel to the axis. Filtered and softened tap water 


was used as the cooling fluid. The inlet and outlet cooling 
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Figure 14. Impingement Cooling Schematic 
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water temperatures were measured along with the cooling 
water flow rate. The heat pipe was rotated in two bearings 
by a three phase variable speed motor with a V-belt drive. 

A small drop-oiler was installed over the main 
bearing housing to lubricate the main bearing. This unit 
consisted of a small reservoir with a needle valve for flow 
regulation. A sight glass permitted observation of the oil 
flow from the reservoir to the main bearing. A flow rate 
of one drop every 20-30 seconds was found to be adequate. 
The end of the evaporator opposite the condenser contained 
a viewport consisting of two glass disks with an air space 
between them to prevent fogging. A strobe light shined 
through this viewport and, when synchronized to the rota- 
tional speed, allowed observation of the boiling action in 
the fluid annulus and the condensate in the condenser. 

2. Instrumentation 

The condenser wall was grooved out externally in 
nine places in Figures 15 and 16. The grooves were spaced 
40 degrees apart and had a depth of 0.010 inch. Welded 
thermocouple beads were soldered into the grooves with a 
Standard lead-tin low-temperature solder. The leads were 
secured to the condenser by banding with twisted wire and 
were passed through the small condenser flange to the mercury 
slip-ring unit on the test stand. Two additional thermo- 
couples sheathed in 0.0625 inch stainless steel tubing 


penetrated the main flange. The sheathes protruded 2.0 inches 
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into the evaporator to measure the vapor temperature and 
were soldered closed at the tip to maintain the system 
vacuum. The inlet and outlet temperatures of the cooling 
water were monitored by Kapton insulated copper-constantan 
thermocouples. All thermocouples were read directly on a 
Newport Model 267 digital pyrometer. Calibration was 
performed in a Rosemount constant temperature bath using a 
platinum resistance thermometer as a standard. A graph of 
thermocouple error as a function of thermocouple reading 
was plotted to be used to correct each experimental reading. 
Cooling water flow rate was measured with a standard rotameter. 
Accurate regulation of the flow rate was accomplished by 
adjusting the pressure regulator installed in the cooling 


line. 


Be EAPERIMENTAL PROCEDURES 
acl asee tone: Condenser Wall 
Tucker [20] established the procedure for prepara- 
tion of the copper condenser for obtaining film condensation 
with water. This procedure was followed with the modifica- 
tions listed below: 
a) The heat pipe was tilted so that the evaporator 
was depressed slightly. 
b) The standard cleaning steps were followed by 
scrubbing with a bristle brush. By depressing 
the evaporator end, cleaning fluids were prevented 


from running back into the heater assembly. 


oe 





c) Complete wetting of the condenser wall was 
observed after cleaning. 
22 ELeorecdieeentang Procedure 

Instead of using the filling procedure of previous 
Studies, a Simpler procedure was used. The condenser was 
depressed, to nearly the vertical, with the evaporator end 
raised up. The heat pipe was then filled with about 350 ml 
of working fluid. The O-ring and glass viewports were 
installed. The bolts of the retainer ring were torqued 
to thirty inch-pounds in ten inch=-pound steps to avoid 
cracking the glass window. This filling procedure was 
followed by a venting procedure, used by Tucker [20], with 
the modifications listed below: 

a) Tilt the heat pipe with the evaporator end 
down, until the working fluid covers up the 
viewport. 

b) Rotate the heat pipe by hand so that the vent 
screw 1S on top. 

c) Slowly apply the power to the heater to about 
0.25 KW until the saturation temperature in 
the heat pipe is at least 220 degrees Fahrenheit 
for the distilled water. 

d) Back out the vent screw allowing some working 
fluid vapor to escape. Control the flow of the 
vapor with the vent screw, being careful to keep 
the pressure in the heat pipe above atmospheric 


pressure. 
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e) Vent for a minimum of ten minutes for distilled 
water, to allow escaping steam to drive off 
alr trapped within the system, leaving about 
250 ml of working fluid left over. 

ff) Turn off the electrical power and reseat the 
VenMe SCEeW . 

g) Return the heat pipe to the horizontal position, 
and turn on the cooling water to the condenser. 

3. Run Procedure 

With filling and venting completed, the system was 

ready to be run. The standard run procedure was as follows: 

a) Open cooling water supply to the condenser and 
bearing cooling water lines. 

b) Start the cooling water supply pump, and adjust 
the rotameter to the desired flow rate. In this 
experiment, only the 50 percent flow rate was used. 

c) O11 the evaporator bearing. 

d) Energize and start the drive motor. Increase the 
speed to 1100 RPM to insure an annulus has formed 
in the evaporator. 

e) Set rotational speed to the desired value. 

£) Let the system reach steady state; take zero 
power thermocouple, cooling water flow rate and 
tachometer readings. 

g) Set power to heater to desired level and allow 


the system to come to steady state. Take 


E> 





tachometer, cooling water flow rate, thermo- 
couple, voltage and ampere readings. 
h) Repeat step g) for each power level desired. 
i) Repeat step f). 
3) Repeat steps e) through i) for each rotational 
speed desired. 
4. Data Reduction 
Data were reduced by correcting the recorded thermo- 
couple readings to obtain a true reading in degrees Celsius. 
Correction was accomplished using thermocouple calibration 
curves. The Celsius readings were then converted to degrees 
Fahrenheit. Percentage flow rate was converted to pounds 
per hour uSing the rotameter calibration curve. Heat trans- 
fer for each run was calculated as the product of the cooling 
water mass flow rate, the difference in temperature between 
cooling water outlet and inlet and the specific heat of 


the cooling water uSing the basic equation: 


QO = mC ,AT (40) 
P 
where 
m = cooling water mass flow rate, lbm/hr 
C = specific heat of the cooling water, taken 
P as one BTU/lbm-deg F 
AT = the temperature difference between the 


cooling water outlet and inlet, deg F. 
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To correct for heat generated by friction in the 
Garlock seals and by viscous dissipation of the cooling 
water, the average zero power heat transfer rate was sub- 
tracted from the measured heat transfer rate in order to 


give a corrected heat transfer rate. 


TET 
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IV. DISCUSSION OF EXPERIMENTAL RESULTS 


A. EXPERIMENTAL RESULTS 

The various experimental runs are shown in Table ll. 

All the heat pipes in this thesis were run at 700, 1400 and 
2800 RPM, except the 1.00 inch diameter cylindrical condenser 
which was run at 1800 and 2800 RPM for the fourth run and 
1400 and 2800 RPM for the fifth run due to erratic readings 
of most of the thermocouples at the low rotation speed. 

It was concluded that the problem existed in the slip-ring 
unit and was due to dirty mercury in the slip-ring. So 

after the fifth run, the slip-ring unit was cleaned by 

taking the mercury and using methanol to clean the inside. 
Nitrogen was blown into the unit to dry it out. New mercury 
was added prior to the sixth run. The plots of heat transfer 
rate versus saturation temperature are referred to as the 
performance curves. 

Figures 17, 18, and 19 show the results of different 
condensers. In all cases, the heat transfer rate increased 
with increasing RPM. In the cases of the rotating cylindri- 
cal condenser, the condensate flow is induced not by centri- 
fugal force but by a hydrostatic pressure gradient established 
in the condensate as a result of its variable film thickness. 
By increasing the diameter of the cylinder from 1.00 inch 
to 1.46 inches, the heat transfer capacity increases. For 


the 0.5 degree truncated cone, the heat transfer is higher 
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Run 
No. 


TABLE ll 


Summary of Experimental Runs 


Geometry 


Gwlaneawacal 1. 
inch diameter 


Cyne tacalk: 14 
inch diameter 


Cyanine eal. 5 
inch diameter 


Cyl TniGmiceal. 1. 
inch diameter 


Gye tid rea); 
inch diameter 


46 


46 


46 


00 


00 


0.5 degree trun- 


cated cone 


0.5 degree trun- 


cated cone 


fe, 


RPM 


700 
1400 
2800 


700 
1400 
2800 


700 
1400 
2800 


1800 
2800 


1400 
2800 


700 
1400 
25010 


700 
1400 
2800 


Test Fluid 


Water 


Freon ers 


Ethanol 


Water 


Water 


Water 


Water 
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Figure 17. Heat Transfer Rate Versus Saturation Temperature 
fox Cylindrical 1.00 Inch Diameter Condenser 
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—-— Data of Wagenseil [19] at 700 REM 


——— Data of Wagenseil [19] at 2800 RPM 


Typical error limits (Appendix B) 


16 


14 


O 700 REM 
4 1400 REM 
Q 2800 REM 





20200) 110° 120 130 140 150 160 170 180 190 200 210 220 230 
Saturation Temperature, T. deg F 


Figure 18. Heat Transfer Rate Versus Saturation Temperature 
for Cylindrical 1.46 Inch Diameter Condenser 
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than the cylindrical condenser. This improvement is due 
to the centrifugal force which increases the axial force 
component on the condensate and thereby increases the con- 
densate flow rate. This increase in working fluid circula- 
tion simply means that more vapor can be condensed per unit 
time. The increase in centrifugal force also flattens the 
condensate film and thereby reduces the overall thermal 
resistance in the condenser. 

The wall temperature profiles as shown in Figures 20, 
21 and 22 were taken from the 1.46 inch diameter cylin- 
drical condenser at different RPM. Figure 23 was taken 
from the 0.5 degree truncated cone at 700 RPM. Distance 
zero corresponds to the condenser end nearest the evaporator. 
The temperatures would be expected to be nearly constant 
if the heat pipe were truly isothermal. For the cylindri- 
cal condenser, the highest temperature is at the evaporator 
end, and the lowest temperature is at about 0.25 inch from 
the evaporator end. For the truncated cone, the tempera- 
ture profile decreases almost linearly, and the highest 
temperature is at the evaporator end. This increase in 
temperature at the evaporator is probably due to conduction 


effects from the heater and the main bearing. 


Pe aNAGYSIS OFTHE RESULTS 
To compare the theoretical and experimental results, 
a log-log plot of the mean Nusselt number versus the Sherwood 


number was made. To make this comparison of theoretical 
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Vem@sUs CxPertimental results, the following simplifications 
were made: 

1) The thermal conductivity of the wall was constant 
along the condenser, 

2) the inside wall temperature was found by using a 
direct application of Fourier's law, 

3) all fluid properties were evaluated at the 
saturation temperature, 

4) heat flux through the condenser was considered 
Lnee Oren, 

5) for the truncated cone, the surface area used the 
average radius between the entrance and the exit 
radius, and 

Gee of OR bh = Q/AAT, AT was taken as the difference 
between the saturation temperature and the condenser 
inner wall temperature. 

The mean Nusselt number for a given heat transfer rate 


1s given by: 


(Q/A AT)L 


% (41) 


Nu = bh Eye 


where 


Nu, = mean Nusselt number 
D = mean heat transfer coefficient, BTU/hr-f£t~-deg F 
A = the average surface area, ae 
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a 
Il 


thetmal conductivity of the fluid, 
BTU/hr-ft-deg F 


L = condenser length, ft 


The mean Nusselt number was determined from the data, 
and this value was plotted against the corresponding Sherwood 
number. Figure 24 shows Nu versus Sherwood number for 
the cylindrical condensers. For water with the 1.46 inch 
diameter cylindrical condenser, the experimental data 
agree with equation (11). The data using the 1.00 inch 
diameter cylindrical condenser, however, show that the Nu, 
increases with decreasing Sh. This trend is the opposite 


to what was expected. The reason for this trend may be 


due to: (1) ripples in the condensate film due to higher 
vapor velocities or (2) a high uncertainty in the low 
power runs. For Freon 113 and ehtanol, with the 1.46 inch 


diameter cylindrical condenser, the experimental results 
are higher than but are parallel to equation (11). The 
reason for this discrepancy is not clear at this time. It 
may be due to differences in vapor densities between water, 
Freon 113 and ethanol, leading to differences in non- 
condensible gas migration, or to differences in vapor 
velocity effects within the theoretical model. 

Figure 25 shows the mean Nusselt number versus the 
Sherwood number of the truncated cone to find out the best 
correlation between the theoretical and experimental results. 


The results of Dhir and Lienhard [8], and Roetzel [14] were 
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chosen to be compared with the experimental results. For 
the one degree truncated cone, using Tucker's data [20] 

and Wagenseil's data [19], the experimental results for 
water fall between both analyses, whereas for ethanol and 
Freon 113, the experimental results are higher than both 
analyses. This may be due again to differences in vapor 
densities between water, Freon 113 and ethanol, leading to 
differences in non-condensible gas migration, or to differ- 
ences in vapor velocity effects within the theoretical 
model. For the 0.5 degree truncated cone with water as the 
working fluid, the experimental results fall lower than both 
analyses, but fall higher than the analysis of Nimmo and 
Leppert Wietersa Straight, cylinder. 

Figure 26 shows Nu_/Nu,, plotted versus the parameters 
SemeectZel seanadiysts. By using all experimental results 
taken to data at the Naval Postgraduate School, it appears 
that good agreement exists between Roetzel's analysis and 
experimental results. The differences of experimental 
results and Roetzel's for the cylindrical condensers are 
less than 15 percent. For the truncated cone condensers, 
the differences are less than 20 percent for water. For the 
truncated cone with ethanol and Freon 113, experimental 
results are higher than Roetzel's analysis by about 35 
percent. However, the agreement is reasonable and the 
analysis of Roetzel can be used to predict rotating heat 
Pipe performance to within + 20 percent for water. For Freon 


113 and ethanol the data lies above the theory by +35 percent. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


mm CONCLUSIONS 


For an internally finned condenser, the analysis of 


the condenser heat transfer problem using a two-dimensional 


wall conduction model and a parabolic temperature distribution 


along the fin can be approximated by the two-element Finite 


Element solution with twenty-five axial increments. 


For a smooth condenser: 


ites 


The theoretical work of Leppert and Nimmo [6,7] can 
be used to predict satisfactorily cylindrical heat 
pipe performance. 

The theoretical equation which best predicts overall 


heat pipe.performance is due to Roetzel [14]: 


23716 237680,,4/ 22 


Sh 


Nu, = Nu [1 + 1.924 tan ¢ ] 


Leo) 
The above equation agrees with the experimental data 
of water to within + 20 percent. The experimental 
data of Freon 113 and ethanol lie approximately 35 


percent above the theory. 


B. RECOMMENDATIONS 


Continue to study the internally finned, rotating heat 


pipe by: 
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Re-examining the Finite Element analysis using the 
mirror geometry presented in Figure 7. 

Analyzing film condensation on a variety of fin 
shapes, and applying the solution of the two- 
dimensional model to determine the shape permitting 
optimum heat transfer rate. 

Manufacturing an internally finned condenser to 
obtain experimental data to compare with the 


theoretical analysis. 


Continue to study the smooth, rotating heat pipe by: 


a 


Using a condenser section with a half cone angle 
greater than 1 degree to obtain additional experi- 
mental data to compare to Roetzel's analysis. 
Examining reasons why the Freon 113 and ethanol 


data are falling above the theory. 
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APPENDIX A 
CALIBRATION OF MEASUREMENT DEVICES 
ine Caletpaaeclon or the Thermocouples 

The thermocouples for both the cylindrical and the 
truncated cone were calibrated, after they had been soldered 
to the outside condenser wall, in a ROSEMOUNT variable 
temperature oil bath utilizing a platinum resistance 
thermometer as the standard of measurement. They were 
calibrated in five degree increments, from 10 to 90 degrees 
Celsius. Each thermocouple was wired to a channel on the 
Newport Model 267 Digital Pyrometer giving individual 
readings in degrees Celsius. For each thermocouple the 


difference T. - T. was calculated where: 


et 
Dae reading in degrees Celsius of the standard, 
and 
_ : : Sela 
T; = reading in degrees Celsius for the 1 


thermocouple at the same temperature. 


The mean of these differences was determined and 
used as the data point for plotting, Ee = De aeeaeas 
versus thermocouple reading in degrees Celsius at known 
temperature. The maximum difference between the mean 
reading and each individual reading was one degree Celsius. 

2. Calibration of the Rotameter 


The cooling water flow rotameter was calibrated 


using a Toledo scale, a large container and an electric 
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timer. Calibration proceeded in 5 percent flow increments 
from 15 to 90 percent flow rate. Water was collected for 
as long as 300 seconds at the lower flow rates and for a 
minimum of 140 seconds at the maximum flow rate. Each flow 
rate was measured two times. The flow rate (lbm/hr) was 
plotted versus percentage flow through the rotameter. This 


curve was used to determine mass flow rate. 
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APPENDIX B 


UNCERTAINTY ANALYSIS 


The method of Kline and McClintock [21] was used to 
estimate the uncertainties in the experimental heat transfer 


rate. The heat transfer of the heat pipe is given by: 


EOra!: 


Qerictional loss Qe le ihe os 
Qcorrected = ae 7 Qe 
where 
Q. = total heat transferred to cooling water, 
EMM OA alg 
Qe = heat transferred to cooling water due to 
trceuon, SrU/AYr 
Q. = corrected heat transfer rate, BTU/hr 
m = mass flow rate of cooling water, lbm/hr 
- = specific heat of cooling water, BTU/lbm-deg F 
AT, = difference of cooling water outlet and 
inlet temperatures, deg F 
AT. = difference of cooling water outlet and 
inlet temperatures due to heat of friction, 
deg F 


ops: 





The uncertainties of these quantities are designated: 


and W 


u AT 


€ 
p 


Wa! 


Wo! 


The fractional uncertainties are given by 


W 








W 
Q W- C W 
se = A + BP 
IG m p 3 
W W W 
Q We c AT 5 WZ 
ee ee 
12 m p 2 
The uncertainty for the corrected heat transfer rate 
given by: 
Wo = (Wo - ee 
C c fs 


The values and 


were taken to be: 


os = 1.000 BTU/lbm-deg F 
Wa = + 0.002 BTU/lbm-deg F 

Pp 
m = 50 percent of maximum flow rate = 
W: = + 1 percent 

m 
W , W = 

ee OS edeq 

= 2 2 Ly as 
Wap = CW + Na ) 
1b} 9 out 
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(Ba) 


is 


U5 ) 


uncertainties of the measured quantities 
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Using the 1.460 inch diameter cylindrical condenser run 


at 700 RPM, the uncertainty limits are: 


T. (deg £) Q. (BTU/hr) UNCERTAINTY, Wa. 
TOR 1611.0 GIG + 39.4% 
133.8 322200 + 637.7 + 19.8% 
LSE 2 222250 + 641.4 + 19.9% 
168.4 5155.2 + 643.6 Sea 
1 Se. 7088.4 + 651.8 + 9.2% 
218.9 T7228 + 655.0 + 3.5% 
2 2idet B87 7.22 + 658.4 ees 
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